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Abstract 
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. Nuclear medium effects in the neutrino cooline; of neutron stars through 

m ■ 

^0 I the reaction channel 77 — > vr*^ — > are incorporated. Throughout 

g ■ the paper we discuss different possibilities of right-handed neutrinos, massive 

fH I left-handed neutrinos and standard massless left-handed neutrinos (reaction 

is then allowed only with medium modified vertices). It is demonstrated that 
multi-particle effects suppress the rate of this reaction channel in the dense 



c/3 

' hadron matter by 6 — 7 orders of magnitude that does not allow to decrease 

- ■ , . existing experimental upper limit on the corresponding tt^uu coupling. Other 

X 

' possibilities of the manifestation of the given reaction channel in different 

physical situations, e.g. in the quark color superconducting cores of the most 
massive neutron stars, are also discussed. We demonstrate that in the color- 
flavor- locked superconducting phase for temperatures T < (0.1 10) MeV 
(depending on the effective pion mass and the decay width) the process is 
feasibly the most efficient neutrino cooling process, although the absolute 
value of the reaction rate is rather small. 



1 Introduction 



Many years ago Pontecorvo and Chiu and Morrison[l] suggested that the process 
77 vv might play an important role as a mechanism for stellar cooling. Gell- 
Mann[2] subsequently showed that this process is forbidden in a local (V-A) theory. 
However, it can occur at the one-loop level which has been computed by Levine[3] 
for an intermediate-boson (V-A) theory, and the stellar energy loss rate through 
77 — >■ vv was found to be smaller than the rates for competing processes (pair 
annihilation e^e' — > vV and photo- neutrino production 76 — > evD). This result 
is not modified when the cross section of the above process is computed in the 
standard model, as was shown by Dicus[4]. Only for very peculiar neutrino coupling 
to photons or unnaturally large neutrino masses this reaction overwhelms the result 
of the standard model[5]. 

There is still another possibility proposed by Fischbach et al. [6], where the 
reaction 77 — > vD could be significant. This is the case when the process is mediated 
by a pseudoscalar resonance and the latter decays into vD due to the existence 
of right-handed neutrinos or due to new interactions beyond the standard model. 
It was assumed that in astrophysical conditions only the pion resonance could be 
important (next in mass not strange ?7-resonance is too heavy in standard conditions) 
and the process was termed the pion-pole mechanism. Thus the process which we 
will continue to discuss in this paper is 77 — > tt*^ — > uV. 

Of course, if the temperature is high enough, and on the other hand, the pion 
dispersion relation in matter allows for the quasiparticle spectrum branch, there 
appears a significant number of thermally equilibrated pion quasiparticles. Then 
the process 7r° — > vD may also be important. In this process the initial thermally 
equilibrated pion is on its mass-shell modified in the matter. In the process 77 — > 
7r° vv the initial reaction states contain no pion, the virtual pion only transfers 
the interaction from thermally equilibrated photons in the 77 annihilation process 
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to produced vv. As we will see below, the process 77 ^ tt" — > has an output of 
the energy Q varying with the temperature as Q oc T** where the power n changes 
with the temperature typically from n = 3 for rather high temperature (T is still 
much smaller than the pion mass = 140 MeV) to n = 11 for low temperatures. 
The process 7r° vv having essentially larger phase space volume (one particle in 
the initial state) yields however the exponentially suppressed output of the energy, 
Q oc 7^3/2g-m,r/T at T < txIt: siuce the initial particle is massive, in difference with 
photons. Concentrating on the discussion of the rate of the 77 — > 7r° — > vv reaction 
channel we shall also compare it with the rate of the competing 7r° — > vv process 
and other relevant processes. 

Temperatures of the order of 10 — 60 MeV are expected in interiors of proto- 
neutron stars formed in supernova (SN) explosions during the early cooling phase of 
the proto-neutron star evolution, according to existing numerical simulations [7]. It is 
exactly in such a situation the pion-pole mechanism is the most effective (T ~ m,r/a, 
a ~ 2 50). Indeed, in Ref.[8] the mechanism was applied to the case of supernova 
SN1987A and it was found to be a quite important process even if the pion partial 
production rate of neutrinos (r(7r'^ ^Bphiyh^R ))) is many orders of magnitude 
smaller than the presently accepted value of the experimental upper limit. In case of 
i^Li^R with the vacuum vertex the process rate is proportional to the squared neutrino 
mass. The result [8] was criticized by Raffelt and Seckel [9] basing on the fact that 
in the calculation of Ref.[8] it was used the vacuum value for the total pion width, 
pmc^ which is a tiny quantity, ~ 0.58 • 10~^m^. They argued that in supernova 
cores pion states are damped mostly by nucleon absorption rather than by the free 
decay so that the width in the medium (F^*^*^) is much larger than in vacuum and 
the pion-pole mechanism should be strongly suppressed. Besides, in the medium the 
process may go on massless left-handed neutrinos through the intermediate nucleon 
particle — nucleon hole states, i.e. as 77 — > 7r° — > nn~^ — > vi,vji reaction channel. 
The processes Tr^j^.^ — > nrT^ — > vjjiji on virtual pions have been discussed in [10]. 
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Without any doubt the neutrino emission from the dense hadronic component 
in neutron stars is the subject of strong modifications due to collective effects in the 
nuclear matter [11], and it is interesting to quantitatively know how much the pion- 
pole mechanism is influenced by the properties of the dense medium and even, if 
the process is strongly suppressed, what is the main cfTcct causing this suppression. 
In this work we compute the cooling rate due to the neutrino emission through the 
pion-pole mechanism in the reaction 77 ^ 7r° ^ i^ri^l (or Uhi^R) for the conditions 
of proto-neutron stars including the effects of dense medium in the pion polarization 
operator (Section 2). Then we discuss other relevant reaction channels. In order 
to clarify their connection to the process 77 — tt'^ ^ i^ri^l (oi' ^l^r) we review 
in the Appendix the optical theorem formalism, see the discussion in [11], that 
allows to calculate consistently the reaction rates including complicated medium 
effects. Continuing to study possible physical situations where the pion pole cooling 
mechanism could be important we consider in Section 3 a possible consequence of the 
pion pole mechanism in the case of neutron stars with the color superconducting 
quark cores. We argue that in the color-flavor-locked superconducting phase for 
temperatures T < (0.1 10) MeV the process is feasibly the most efficient neutrino 
coohng process, although the absolute value of the reaction rate is rather small. 
Then we draw our conclusion. 

2 Emissivity of the process 77 ^ 7r° ^ vv from 
the hadronic matter 

The cross section of the process 77 — > 7r° — > I'Ri'Li^L^R) in vacuum is given by 
.ac(. _ 87r(5Vm^)r(7r° ^ 77)r(7r° ^ v9)F{s) 

" ^ ^ {s- miy + (m,r--)2 ' ^ ' 

where F(s) is an unknown function of s =(c.m. energy)^, representing the product 
of the vertex function for the off-shell processes 77 — > 7r° and 7r° —>■ uu, constrained 
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to F{s = ml) = 1. Following Ref.[6] we may assume F{s) ~ 1 off mass-shell, 
that seems to be a reasonable approximation for the energies and momenta that 
we are dealing with. r(7r° — > 77) is the partial width of the pion decay into 
two photons. Here we assume that r(7r° 77) is approximately equal to the 
total pion width (rjj."'^), for which we use the experimental value. The pion par- 
tial width into neutrinos r(7r'^ — > uu) has the following experimental upper limit 
r(7r° vv)/V'"^°^ < 8.3 X 10~^ [12]. In (1) one recognizes the free pion propagator 
modulus squared | D^""^ p entering squared matrix element of the reaction under 
consideration. 

Strictly speaking, in application to the nuclear medium all the terms in (1) should 
be modified. E.g., the partial width r(7r° 77) may increase with the tempera- 
ture [13]. The pion partial width into neutrinos r(7r° vD) is determined by the 
corresponding 7r° — vD coupling. The squared vertex is averaged over the neutrino- 
antineutrino phase-space volume, see Appendix. If we knew the explicit form of the 
coupling we could explicitly calculate r(7r'^ — > vv) and its temperature dependence. 
However the main modification comes from the change of the pion propagator in 
dense nuclear medium due to the pion pole. Therefore, below we consider only this 
modification. Then, one should replace in (1) the free pion propagator D'^"^ by the 
in-medium one 

^ [a;2_^2_^2_n^„(^,A;,p,y,T)]' 
where Ti^o{w,k, p,Y,T) is the total retarded polarization operator of the neutral 
pion, dependent on the pion energy a;, momentum fc, baryon (nucleon in our case) 
density p, isotopic composition Y = Z/{N + Z), and the temperature T. 

If photons are in thermal equilibrium, the energy loss rate of the process 77 —>■ 
7r° — > i^rVl^ulVji) occurring in dense interior of the proto-neutron star, at which the 
energy is converted into neutrino pairs, is given by 

^ - (2^ J [exp{u:,/ksT) - 1] [exp{u:,/ksT) - 1] + ^^^"-'^'^ ' 
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where uji, lj2 are the photon energies, ki and k2 are their momenta, 9 is the angle 
between photons, Vrei is the relative-velocity factor 

Vrel = \{kl ■ k2 - U;iUJ2) /'^l'^2\ = 1 - COS^, (4) 

and a'^^^ is the medium dependent cross section of the process given by (1) with 
jjvac i-gpiaced by D^o- Eq. (3) is the straightforward generahzation of the result 
[6, 8] obtained with the help of the replacement o"^"'^ ^med ^ 

The energy loss rate (3) can be presented in the following form 



T(7r° ^77)\ /r(7r' 







vv] 



Qmed _ 4 X 10- ^ T^^Ut), erg/icm' ■ s), (5) 

\ J \ m-j, J 

= 140 MeV, Tg — T/IO^K, r — {kBT/rriT^) and I{t) is given by the dimensionless 

integral 

^'^^ ~ Jo e^i - 1 io _ 1 ^ ^ 

'■2 du{xi + X2)u^ 



X 



/ 

J K 



{2xiX2UT^ - 1 - m-meU^of + (m-ir^)^' 
where k = -^^{my/kBTY, rUi, is the neutrino mass, xi^2 — '^i.ijkBT ^ u — \ —cos9, 
m~^T.„ = m~^T^°''^ —m~'^ImIl^o- Further we assume that rriu/kBT <^ 1 and then put 
the lower limit (k) of the integral in the variable u equal to zero, that corresponds 
to using nil, = in all the phase space calculations. However we take rrii, ^ 
into account evaluating the r(7r° ^l^b) width. For the further convenience we 
assume that -Rell^o and Imll^o are the real and imaginary parts of Il^o{w, k,p,Y, T) 
describing only the strong interaction processes. Therefore we separated in (6) the 
value rj^"^ which is due to the electromagnetic and the weak interaction. With J(r) 
replaced to P'^{t) (putting -Rell^o and /mll^o to zero) we reproduce the results of 
Refs.[6, 8]. One easily finds the corresponding asymptotic expression 

/ ^ r-(r « 0.1) ~ 23040 C(6)C(5) [l + 96r\(7)/C(5)] , (7) 

C is the Riemann function, C(5) ~ 1.037, ((6) ~ 1.017, C(7) ~ 1.008. In order to get 
(7) we dropped and i^ell^ and expanded (6) in the value 2xiX2Ut'^ <^ 1. Since 



typical values Xi ~ 2:2 ~ 6 in the resulting integral, the limit expression is valid 
for T ^ 0.1. As we checked numerically the limit (7) is actually achieved only at 
T < 0.01, if is as small as F^"''. Using Eq. (7) and the evaluation of r(7r'^ 77) 
we obtain 

Q^'^'ir < 0.01) ~ 0.6 • 10'° (m;'r(7r° ^ i/P)) erg/{cm^ ■ s) (8) 

at small temperatures. 

When the temperature increases the denominator becomes to be near the pole. 
Then dividing and multiplying /(r) by and using the corresponding presentation 
of the 5-function we roughly estimate 

lf(^- V/^>.>n-A 0-6(e-^/('-)|ln(2r)| + l)e-V('-) 

Prom the latter estimate we recognize the resonant character of the rate. I'""-'^{t ~ 
0.1) is ~ 10^ times larger compared with /^"^(r = 0) for ~ V"^. 
For higher temperatures assuming 2xiX2Ut'^ ^ 1 we estimate 

/ (r » {^f") ^ 12C(3)C(4)/r^ C(3) ^ 1-202, C(4) ^ 1.082, (10) 

that produces g™'^ ~ Q'"^'^ ~ 2.6 ■ lO'^ {m-^V{'K^ ^^))Tl, erg /{cm' ■ s). 

In the so called "standard scenario" of the neutron star cooling one considers the 
modified Urea process as the most efficient process. The emissivity of the modified 
Urea process is estimated [14] with the help of the free pion propagator as Q^^ ~ 
lO'^Tg^ erg /{cm? ■ s). Comparing it with (5), (9) we see that for T > 10 MeV the 
process under consideration would be much more efficient process than the modified 
Urea process, if the pion width and the mass in medium were not essentially changed 
compared to the vacuum values. However in reahty the width F^^"'^ is replaced to 
a much larger medium value. Also the pion mass is modified by the polarization 
effect (thereby the pion pole begins to manifest itself at T > rrfj^ ja^ rather than 
at T > rriT^/a). Thus one may expect that with taking into account of the medium 



effects the estimation (7) is not essentially modified whereas the value (9) must be 
significantly suppressed mainly due to the suppression of the width. 

Techniques for the description of collective effects in dense hadronic matter have 
been developed in the last decades[15, 16]. E.g., the medium effects appearing in 
the pion propagator have been discussed in detail in Rcf.[17] (see Appendix B of 
that work). For the temperatures and pion energies with which we are concerned 
the main modification of the pion polarization operator is due to the density de- 
pendence rather than the temperature dependence and it is, thereby, a reasonable 
approximation to put T in n^o(ci;, fc, p, Y, T). Depending on the values of the 
typical pion energy and momentum different terms can be important in the pion 
polarization operator. At small pion energies Lo/m.,^ <^ 1 and for typical momenta 
k ~ Pfn, Pen is the Fermi momentum of the nucleon, the nucleon particle - hole 
contribution is attractive and the dominant one, what results in the softening of the 
virtual pion mode with increase of the density and leads to the possibility of the 
pion condensation at p > pc > Po, cf. [15, 16, 17]. In our case c<j > A;, as follows from 
the reaction kinematics, and the nucleon particle - hole contribution is minor. Then 
the main terms in the pion polarization operator are the A-particle - nucleon hole 
part and the regular part related to more complicated intermediate states. Thus, 
n^o — + n^gg, and the partial contributions are given by [17]: 

Re U^iu, k, p, y, T = 0) ^ . (11) 

in units ^ = c = 1, ~ 2.0 mT^{p/ pQ)V'^^^k'^Lijis,{t), p is expressed in units po, po — 
0.5 is the saturation nuclear density, the form factor F^^^ Ri F^^^//? ~ 
for the rather small momenta of our interest, /3 ~ 1 + 0.23A;^/m^ is an empirical 
factor taking into account a contribution of the high-lying nucleon resonances. The 
nucleon-A isobar correlation factor T{g'^) is given by 



1 + 



c 



(12) 
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with C ~ 0.9 m^(p/po)r^ArA and 

^4*)-2.1,„.(l + ?i|^)+^. t = e-.\ (13) 



m*j^{p) is the effective mass of the nucleon quasiparticle, 



/.nn^K^,p,F,T^O).- ^f£(^^y , (14) 



7o/i;^ takes into account the A-isobar width, with an empirical value 70 ~ 
0.08 T^^^. In the numerical evaluations below we for simplicity assume 
iij\{t) ^ uj'^ that is a reasonable approximation for the conditions we are dealing 
with. 

The regular part of the pion polarization operator is yet more model dependent. 

Its value is recovered with the help of the pionic atom data, cf. [18], and a procedure 
of going off mass-shell. We present it in the following form, cf. [17], 

Re n^e«(^, p,Y,T = 0)c^ (-0.25 + 0.25 ml + 0.5 k'^) — , (15) 

7mn,«,^(a;,A;,p,y,r = 0)~-0.15m2 l^-^j ^ 

0.09(p/po)^Pe .... 

[l + 0.38(p/po)- 0.047(p/po)^0]'' ^ ' 

The factor ^ which we inserted in these expressions compared to those of [17] takes 
into account asymmetry of the isotopic composition of the proto-neutron star mat- 
ter. Near the pion mass-shell, with ^ = 4y(l — Y) we approximately describe 
results given by the potentials I-III used in [18] for the pion atoms and with 
^ = 1, the results for the potential IV. For the value Y ~ 0.4 typical for initial 
stage of proto-neutron star cooling in both mentioned cases one can put ^ ~ 1. 
The real and imaginary terms shown above are the ones which enter Eq.(6), be- 
ing responsible for the density effects, where we put ^ ~ 1. We also used that 
n^oiou, k,p,Y,T) = I [U^.{u;,k, p,Y,T) + n^+{-io,-k, p,Y,T)], and the linear in 
uj terms entering Rell^± for F 7^ 1/2 do not contribute to ReU^o- 



We computed the energy loss rate including medium effects in the pion polar- 
ization operator. In Figure 1 we present the results for the ratio of the energy 
loss rate (Q™^'^) due to the pion-pole mechanism calculated with medium effects 
taken into account in the pion polarization operator to the one (Q^"'^) computed 
with ReU^ = and the vacuum pion decay width (cf. (8)), for three values of 
the nuclear matter density: p = (0.5; 1; 2)po- We used ~ (0.9; 0.85; 0.7)m]\f for 
those densities and put m,r — 140 MeV. The ratio Q"^^'^ / Q'""-'^ is depicted in Figure 
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Figure 1: Ratio of energy loss Q'^'^'^ / Q^°-'^ as function of temperature for p = 0.5po 
(solid line), p = Po (dot-dashed line) and p = 2po (long dashed hne). 

1 as a function of temperature. This ratio does not depend on an unknown value 
of r(7r° up). We see that the nuclear matter effects decrease the output of the 
energy typically by six to seven orders of magnitude depending on the temperature 
and the density, in agreement with the expectations of Raffelt and Seckel [9]. We 
also show in Table 1 the results for the dimensionless integral of Eq.(6). 

Qualitatively, the reasons for the enormous decrease of the energy loss rate are 
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Table 1: Numerical results for the dimensionless integral /(r) of Eq.(6). 

related to the strong pion absorption in nuclear matter, as has been predicted in 
Ref. [9]. The total pion width at appropriate pion energies and momenta grows 
up to tens of MeV with the density, that is orders of magnitude larger than the 
vacuum contribution F^""^. Other reason is the energy-momentum dependence of 
the real part of the pion polarization operator. This is also not a small change 
because the mechanism is totally dependent on the resonant behavior, and outside 
the resonance the contribution sharply decreases. We do not need to be on the 
top of the resonance, however we cannot be too far away either. Medium effects 
push the typical pion energy to a larger value {Rell^{u!, k) > at pion energies and 
momenta in the region of the pole) . 

The ratio r(7r° z/z/)/r(7r° 77) in many models of elementary particles 
beyond the standard one is proportional to the neutrino mass {mu). In Ref. [8] a 
limit on the neutrino mass was obtained due to the strong constraint on r(7r° uP). 
However, with medium effects included into consideration no astrophysical limit 
on m,^ better than the existent ones can be obtained. A nonzero neutrino mass 
induces other mechanisms much more efficient than the one generated by the pion 
resonance in this situation, and these mechanisms provide then tighter constraints 
on the neutrino masses (see, for instance, Ref. [19]). 

Above we discussed only the role of the reaction channel 77 — > 7r° — > i/P. But 
there are other competing reaction channels. The relation between all these pro- 
cesses becomes to be clear if one apphes the so called optical theorem formahsm for 
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the calculation of the reaction rates, see discussion in the Appendix. 

First, if the pion can be described within the quasiparticle approximation in 
some region of its energy and momentum, one has a contribution of the process 
7r° vV whose rate is given by, 

Q -[-^^ )n.m^, (17) 

Here we used the dispersion relation u"^ — (m^.y + vlk'^ for small momenta k ~ ksT 
typical in this reaction. The effective mass m* ~ and the velocity < 1 are 
calculated according to eqs. (11) - (15). The effects of the partial r(7r° —>■ 77) 
width are not present in this process. 

Second, the presence of finite width of the virtual pion /mll^o ^ due to the 
strong interaction means also a finite contribution of the reaction 7r°j^j — > vP that 
does not need 77 and pion quasiparticle states, but relates to the corresponding nu- 
cleon states, cf. [10, 11] (it is clearly seen after the cut of the in-medium pion Green 
function describing propagation of the in-medium pion) . Presence of the imaginary 
part of the nucleon-hole term of the pion propagator (with the full vertex) in ap- 
propriate energy-momentum region would lead to a contribution of the processes 
N Nuu, NN NNuu going via the virtual pion, the later couples with 
pu. The presence of the imaginary part of the regular term in the pion polarization 
operator corresponds to more involved multi-nucleon and multi-pion states. 

Third, the processes as 77 — > 7r° — > nn~^ — > i-'iP^ {n~^ is the neutron hole) with 
ordinary left-handed massless neutrinos are possible leading to substantially larger 
contribution to the emissivity than that related to the massive left-handed neutrinos 
in the reaction 77 — >■ 7r° — >■ i'li^r. The reaction channel 77 — >■ 7r° — >■ nn^^ i^l^r 
may also lead to substantially larger contribution to the emissivity than that related 
to the right-handed neutrinos (depending on the value of the r(7r'^ — > i^rVl)) in 
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the reaction 77 — > tt" — > z/rz/^ which we have considered. Besides, the process 
7r° — > nn~^ — > i^Lp^ going on the left-handed massless neutrinos is possible going on 
the thermal pion quasiparticle. The latter rate can be calculated in complete analogy 
to that computed for the massive pseudo-Goldstone (photon) mode in [20]. The 
processes involving pion but going via nn^^ — vv coupling (with usual V —A coupling 
of N to vlt^r) exist only due to the medium modification of the vertices. Another 
possibility, which will be discussed below, is that due to the breaking of the Lorentz 
invariance in matter there might appear two pion decay coupling constants, the 
temporal one, /t, and the space one, /s,[21, 22]. Their finite difference would also 
lead to a contribution to the 7r° coupling with the left-handed neutrinos. Although 
all these relevant processes are related to each other, they are characterized by 
quite different phase-space volumes, kinematics and/or vertices. Having selected 
the 77 — > TT*^ — > channel among others, we considered the enhancement of the 
rate due to the particular kinematics of this process (massless particles in initial 
states). But, as we have shown, the presence of not as small pion width in the 
hadron matter significantly suppresses the rate. 

Concluding this section we once again stress that although the temperature 
dependence of the energy output in the process 77 — > 7r° ^ uV is quite different 
from those for the other processes, where also enters the t^vV coupling, the absolute 
value of the rate is proved to be strongly suppressed due to the presence of not as 
small pion width in the dense and heated nucleon matter. Thus, due to the width 
effects, the process has a minor role in the cooling of the dense hadron matter. 
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3 Emissivity of the process 77 ^ vu from 

the color-flavor-locked phase 

There is another interesting possibihty in connection with the process under consid- 
eration. As has been recently shown, the interiors of the most massive neutron stars 
may contain dense quark cores which are high temperature color superconductors 
with the critical temperature Tc ~ 0.6Ag < 50 MeV, where is the pairing gap 
between colored quarks, see the review papers [23, 24]. Also the possibihty of self- 
bounded strange quark stars, being the diquark condensates, is not excluded, see 
[25]. The diquark condensates may exist in different phases. The most symmetric 
phase of dense quark matter is the so called the color-flavor-locked phase. This 
phase becomes to be energetically preferable in the large density limit. The neu- 
trino processes are significantly suppressed in this phase due to presence of the large 
diquark gap and absence of the electrons [26, 27]. On the other hand, it was shown, 
cf. [28, 24], that this phase contains low- lying excitations with the pion, kaon and 
77, 77' quantum numbers (m*o, m^o, ?77.^o, "^*, m*/ are in the range (1 — 100) MeV, all 
in-medium masses are indicated by m*). Excitation spectra, as they are calculated 
in the mentioned works, contain no widths effects, the i-meson Green function is 
given by Di = ^^2_fc2/3,^*:i) ■ However, in any case there is a width contribution at 
least from the process tt^ 77. Then —ImU^ ~ r(77 — >■ tt^) ~ r''"^(77 — >■ 7r°). 
Although more involved effects may also simulate the corresponding width terms, 
one may expect that the width effects are, nevertheless, rather suppressed due to the 
nature of the superconducting phase with the large gap. In this situation the pion 
pole mechanism could become an efficient cooling mechanism. Also other meson 
poles may essentially contribute. 

For the process under consideration the emissivity of the superconducting 
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medium Q*^*" is given by expressions similar to Eqs. (3) to (5). We take 

' {u^-ky3-mll)^ + {ImURy ^^^^ 

and replace it into (1) instead of (2). Other quantities in (1) are assumed to be 
unchanged. Then the energy loss rate is presented as 



Q^^ c 4.2 X 10- flV^") fn^'^--)] Tl\''-^)H^^{fl 



(19) 



now with the extra factor (^)^ and the integral /(r) replaced by 

xf dxi f°° X2 dX2 



^ ^ io e^i - 1 Jo 



e^2 _ 1 



du{xi + X2)U^ , . 

^ A {lf^[{xi + X2)^ + x,X2u]-iy + {-ml-^ImU^oy' ^ ' 



T — T/m*^. Again we further take k = 0. For r <^ 0.1 we get 
/^^(r < 0.1) ~ 23040 C(6)C(5) 



l + 92r^ fZCa + OT)' 
+ ^23C(6) + C(5), 



f211 



that not essentially differs from the estimation given by Eq.(7), C(8) = 1.004. 
The value (21) is independent on m* for r 0.1, and Q^^ oc T^^(m*)"^, 
QSCjQvac ^ {m^/mlf, cf. (8). For [ml'^ / {-IQ'^ImU.^)^/^ > f > l/a we 
roughly estimate I^^ ~ 0.1[-m;-2lmn^]-if-^(e-V(2?)|ln(2f)| + l)e-V(2T)^ ^ ^ nfi 
for -Imn^ ~ r^^^ and oc {mlf e"^* ' ^'^'^^ {-luiYi^) . For still higher 

temperatures J'^*" ~ and Q^'~^ oc T'' being almost independent on m*, 

qSC ^ 3 . 2o29 (77j,-ir(7r° ^ i/p)) rj', erg/{cm^ ■ s). Interpolation estimation being 
roughly valid for all temperatures produces 

^sc _ 24299(1 + 120r'e-^°") 0.l(e-V(2^)|ln(2f ) | + l)e-^/('") 

~ 1 + 2880?^ ^ f8[-m;-2Imn«] ' ^ 

The emissivity {Q^'^) computed with the above interpolation formula (Eq.(22)) and 
for the total pion width given by the vacuum one is shown in Fig.2. The value 
r(7r'^ — > up) is assumed to be equal to its experimental upper limit. This figure 
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serves as a guide for the numerical calculation of Eq.(20), showing the very fast 
increase of the emissivity as we approach temperatures near the pion pole one and 
demonstrating a slow increase as we go to temperatures above the effective pion 
mass scale. The interpolation formula and the full numerical calculation of the 
emissivity agree with high accuracy for small T (compared to m*) and differ by a 
factor < 2 for T ^ m* . 

In Figs. 3, 4 we show the numerically calculated emissivity of the process under 
consideration for the superconducting media, Q^'^ , in the temperature range from 
5 to 50 MeV for m* = 10 MeV and m* = 50 MeV, respectively. These curves were 
obtained assuming r(7r° vv) equal to its experimental upper limit. In the given 
temperature interval (from 5 to 50 MeV ) the emissivity curves are scaled up for 
m* =50 MeV compared to m* = 10 MeV in accordance with above estimations 
of I^'^ and interpolation equation (22). The pole asymptotic manifests itself even 
in the case when the width is rather suppressed (solid curves). It happens for 
T < (15 ^ 20) MeV in case m* = 10 MeV and in the whole temperature interval ( 
from 5 to 50 MeV) in case m* = 50 MeV. For m* = 10 MeV solid, dash and dotted 
curves reach the high temperature asymptotic for T > 20 MeV. For m* = 50 MeV 
the solid curve achieves the high temperature asymptotic for T > 50 MeV. The 
smaller is the width the more pronounced is the pole mechanism. This is clearly 
seen if we compare Fig. 2 (the result of the interpolation formula computed for the 
vacuum width) with Figs. 3, 4 where we assumed that more involved effects may 
simulate a width larger than the vacuum one. 

Compared to the hadron case, in the color-flavor-locked phase of the supercon- 
ducting quark medium the effective pion mass is much smaller and the width is 
suppressed. In this case the resonance effect appears with a larger amplitude (due 
to very small width) and at smaller temperatures. Thus we argue that for the star 
core in the color-flavor-locked phase at T < T^^, i.e. beyond the neutrino opac- 
ity regime, in absence of any other efficient cooling reaction channels the pion pole 
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mechanism could become an efficient cooling mechanism even if r(7r° — > i/i?) being 
by many orders of magnitude below the value determined by the modern exper- 
imental upper limit. Here the value of T^^^ is different from that for the usual 
proto-neutron star being determined by the condition that the neutrino mean free 
path in the quark core becomes to be comparable with the size of the core. In 
the neutrino opacity regime, T > T,^^^., the reaction — > delays the neutrino 
transport from the quark core to the hadron shell. To come to these conclusions 
we used only that the effective pion mass in the CFL phase is essentially smaller 
than the vacuum one, the pion width is small and the pion momentum is shifted by 
k A;/V3. In some relevant temperature interval for the width ~ F^"^ and for 
not too small effective pion mass the rate Q^^"' may even exceed the cmissivity of the 
most efficient direct Urea process being Q^^ ~ 10^^(T/10 MeV)^ erg/{cw? ■ s) for 
the non-superfluid (hadron) neutron star matter. Also at very low temperatures the 
emissivity of the process decreases with the temperature according to the power low 
{Q^*" oc T^^) rather than exponentially as other relevant processes in the CFL phase. 
Therefore in spite of a low absolute value of the rate at such low temperatures, the 
process is the dominating process. 

Now we may compare the emissivity of the process 77 — > 7r° — > vV with the 
estimation for the emissivity (17) of the process 7r° vv in the quark matter 
discussed in [26, 21, 22]. The emissivity of the later process is given by Eq.(17) 
with = 1/3 for the color-flavor-locked phase. We see that the ratio (19) to (17) 
is larger than unit in a wide temperature interval of our interest, ksT ^ 7 MeV for 

~ 70 MeV and /csT < 0.2 MeV for m*^ ~ 10 MeV. In this estimation we again 
used that the total pion width is almost exhausted by the 77 decay, F(77 — > 7r°) ~ 
—lmX\^ and, as before, we took F(77 — >• 7r°) ~ F^"^. If the width F^^ were larger 
as the consequence of the presence of some processes which were not considered up 
to now, then the pole contribution could be more suppressed (see dotted, dash and 
solid curves in Figs 3, 4) and the process 77 — > 7r° ^ vv would be relevant only for 
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low temperatures T < (0.1 ^ 1) MeV. 

Notice that the ratio of the reaction rates for 77 — > 7r° — > vV and 7r° — > vV does 
not depend on the value of r(7r° — > uV) and it always becomes larger than unit with 
the decrease of the temperature (however the critical temperature when the rate 
reaches the unit depends on the vahies of the parameters). 

Refs [21, 22] found an interesting possibility that in the color- flavor- locked phase 
the process — > vV is allowed also in the standard model with left-handed neutrinos 
even if the neutrino mass is zero. The important point that was noticed is that the 
temporal and thermal components of the pion decay constant need not be the same 
in the dense medium due to the violation of the Lorentz invariance. The calculation 
of [21] yields 

- °-^-^°"iiow)^ioi7(i(^^^ 

jiq is the quark chemical potential. We see that for m* ~ ksT ~ 10 MeV the 
value (23) is by three-four orders of magnitude smaller than the value of the exper- 
imental upper limit for r(7r'' uP) for the right-handed neutrinos. Therefore both 
possibilities should be studied. 

Ref. [21] also evaluated the value r(77 — > 7r°) for the corresponding decay in 
the color-flavor-locked phase. This estimation only by the factor ~ 1 differs from 
the value r(77 — > 7r°) ~ F^""^ which we used in our estimations. 

4 Conclusion 

In conclusion, (i) we evaluated the cooling rate of the proto-neutron stars via the 
pion pole mechanism taking into account the nuclear medium effects and we showed 
that with the width effects included this mechanism gives no constraints on the 
corresponding 7r° — > ^'r^'l(i^li^r) decay width, (ii) We also discussed possible con- 
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tribution of this mechanism to the coohng of different astrophysical systems, as 
is the case of neutron stars with dense quark cores. We found that the process 
77 — > 7r° — > yV is proved to be the most efficient process in the color-fiavor-locked 
superconducting core at temperatures T < (0.1 -=-10) MeV depending on the effective 
pion mass and the pion decay width. Depending on what the future tells us about 
r(7r° w) and on the possibility of the color superconductivity in neutron stars 
we believe that these questions deserve a further more detailed study (perhaps even 
further studying of a dependence on medium effects of all the terms appearing in 
the cross section of 77 — > tt*^ — > ^'rPl(^'lPr)). 
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A Optical theorem formalism 

In [29] and then in [30], see [11] for a review, there was developed the optical theorem 
formalism in terms of full non-equilibrium Green functions to calculate the reaction 
rates including finite particle widths and other in-medium effects. Applying this 
approach, e.g., to the antineutrino-neutrino production we can express the transition 
probability in a direct reaction in terms of the evolution operator 

dWT_^,, = ^^y^'^ V<0|^t|p^ + x><Pi/ + X|^|0>, (24) 

where we presented explicitly the phase-space volume of Pv states; lepton occupa- 
tions of given spin are put equal to zero for 1/ and P which are supposed to be radiated 
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directly from the system. The bar denotes statistical averaging. The summation 
goes over the complete set of all possible intermediate states {X} constrained by the 
energy-momentum conservation. Making use of the smallness of the weak coupling, 
we expand the evolution operator as S ^ 1 — i Jl^^ T {Vw{x) Snuci{x)}dxo ■, where 
Vy/ is the Hamiltonian of the weak interaction in the interaction representation, 
S'nuci is the part of the S matrix corresponding to the nuclear interaction, and 
is the chronological ordering operator. After substitution into (24) and averaging 
over the arbitrary non-equilibrium state of a nuclear system, there appear chrono- 
logically ordered {G ), anti-chronologically ordered (G"*""*") and disordered (G^" 
and G ^) exact Green functions. 

In graphical form the general expression for the probability of the neutrino and 
anti-neutrino production is as follows 




representing the sum of all closed diagrams (— in~+) containing at least one ( — h) 
exact Green function. The latter quantity is especially important. Various contri- 
butions from {X} can be classified according to the number N of G lines in the 
diagram 



V V V V 




V V V V 



Being expressed in terms of the exact Green functions, each diagram in (26) repre- 
sents a whole class of perturbative diagrams of any order in the interaction strength 
and in the number of loops. This procedure suggested in [29] is actually very help- 
ful especially if the quasiparticle approximation holds for the intermediate fermions 
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or bosons. Then contributions of specific processes contained in a closed dia- 
gram can be made visible by cutting the diagrams over the (H — ), ( — h) lines. In 
the framework of the quasiparticle approximation for the non-relativistic fermions 
Gp^ — 27rinFS{e+iJ,—£p—ReT:^{e+iJ,,'p)) {np are fermionic occupations, for equilib- 
rium Up = 1 / [exp{{e — /jp) /T) + 1]) , and the cut eliminating the energy integral thus 
requires clear physical meaning. In this way one establishes the correspondence be- 
tween closed diagrams and usual Feynman amplitudes although in the general case 
of finite fermion width the cut has only a symbolic meaning. Next advantage is that 
in the quasiparticle approximation any extra Gp^, since it is proportional to n^, 
brings a small {T/sp)^ factor to the emissivity of the process. Dealing with small 
temperatures one can restrict by the diagrams of the lowest order in {Gp^Gp~), 
not forbidden by energy-momentum conservations, putting T = in all G"*""*" and 
G Green functions. For the relativistic bosons with the non-conserving number, 
in the quasiparticle approximation = — 27rms5 (00"^ — m% — k"^ — Ren''^(ct;, /c)) 
(ub are bosonic occupations, for equilibrium ub = 1 / [exp (cj/T) — 1]). In a wide 
temperature and energy-momentum region Imll^ is not small and the quasiparticle 
approximation is then not valid for the boson Green functions whereas the region, 
where ImE^ is small and the quasiparticle approximation is valid for the fermion 
Green functions, is usually much wider. 

If one is interested only in the processes related to the uu coupling with the 
7r°, the hatched block in (25) is reduced to the exact D'o' Green function. This 
Green function satisfies the exact Dyson equation. Only if in a specific region 
of the pion energies and momenta {cu and k) the pion Green function can 
be approximated by the 5-function, integrating over this region one may use the 
quasiparticle approximation for the pion. In such a way one usually calculates 
the rate of the 7r° — > i/P process. Certainly, in other regions of pion energies and 
momenta the pion Green function contains the width relating to different channels 
of the pion decay. E.g. the polarization operator of the 7r° contains the " — h" 
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77 loop. This term corresponds to the contribution D't^ G1^~ D^o' . Within the 
quasiparticle approximation to the 7 one cuts the hues and gets in this 

way the contribution of the 77 — > 7r° — > uV process which we discuss in this paper. 
When considering right-handed neutrinos we do not know the exphcit expression 
for the Ti^ vv vertex. Several different expressions can be used. Therefore we will 
express the result of the integration in the vv states in (24) via the phenomenological 
value of the width r(7r° vv) for which there exists the experimental upper limit. 
If we knew the coupling we could present an explicit calculation, as one usually does 
for the left-handed neutrinos. 
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Figure 2: Emissivity Q'^*^ in erg/{cw?' • s) computed with the interpolation formula 
of Eq.(22) for m* = 10 and 50 MeV. The curves were obtained assuming the total 

width ~ r^"'^ . 
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Figure 3: Emissivity Q*^*^ in erg/{cm^ ■ s) as function of temperature for 
—ImU^/ml ^ 1 MeV (solid line), —ImU^/ml ~ 0.1 MeV (dashed line) and 
-ImU^/ml ~ 0.01 MeV (dotted line) for ml = 10 MeV. 
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Figure 4: The same as in Fig. 4 but for m* = 50 MeV. 



26 



